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O , Abstract 

, A non-parametric statistical model is constructed to directly relate the 

' distribution of observed microlens timescales to that of the mass function 

of the population from which the lenses are drawn, corrected for obser- 
vational selection based on timescales and event amplifications. Explicit 
distributions are derived for microlensing impact parameters and maxi- 
mum amplifications; both are shown to be statistically independent of all 
other parameters in the problem, including lens mass. The model is used 
to demonstrate that the narrow range of microlens timescales observed 
toward the Large Magellanic Cloud (LMC) is probably not consistent 
with lensing by a widely distributed spheroidal population ( "dark halo" ) 
' of large velocity dispersion, but is consistent with lensing within a rotat- 

. ing thick disk. Poor numerical conditioning of the statistical connection 

between lens masses and event timescales, and small number statistics, 
severely limit the mass function information obtainable from current mi- 
crolensing surveys toward the LMC. 



o 

^ : 1 INTRODUCTION 

Microlensing is a promising tool for eliciting the properties of massive compact 
objects in the Galactic halo and toward the Galactic bulge, and several such 
, observing programs are underway or have been completed. In particular, the 

' MACHO and EROS projects have published results for several years of mon- 

itoring > 10^ stars in the Large Magellanic Cloud (LMC; Alcock et al. 2000, 
Milsztajn et al. 2001), which detected 10-20 events; and both the MACHO and 
OGLE projects (Udalski et al. 2000) have published the results of microlensing 
surveys toward the Galactic bulge that include more than 200 events. Several 
additional microlensing surveys are currently underway. The results of these 
surveys have provided otherwise inaccessible information on the structure and 
contents of the Galaxy, but not always in unambiguous ways. The primary 
object of interest has been the total mass in the form of compact objects, par- 
ticularly in the (presumably dark) Galactic halo, usually analyzed in the form of 
microlensing optical depths originally proposed by Paczyhski (1986). As useful 
as these analyses have been, little information on the mass function of the popu- 
lation underlying the microlensing observations has been forthcoming, although 
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some parametric modeling (e.g., Green 2000) has sufRced to indicate some pos- 
sible constraints conditional upon various kinematic models of the halo. Here 
we propose a more direct way to elucidate the mass function of the population 
responsible for microlcnscs, in the form of a non-parametric statistical model 
that directly relates the mass function to the microlens quantities being ob- 
served. In principle, such a model should allow an unambiguous determination 
of the mass function, given only the kinematics of the population. In practice, 
we expect that the information obtainable from such a method will be limited 
by small number statistics and ambiguities concerning the kinematics of the 
lensing population. 

The necessary microlens physical relations are as follows (see, e.g., Paczynski 
1996 or Gould 1996 for more complete discussions). The physical scale of a 
microlensing event in the lens plane (i.e., the plane of the sky at the lens location) 
is given by the Einstein radius: 



2 GMDl {Ds-Dl) 
^-=cV ' 

where and Ds are, respectively, lens and source star distances from the 
observer; and M is the lens mass. For solar-mass halo microlenses viewed against 
the backdrop of the LMC, te is typically a few astronomical units. In terms of 
this quantity, the scaled impact parameter (3 and event timescale T are defined 
to be 

= b/rE. (2) 
T = teIV^, (3) 

where &, the "impact parameter", is the closest approach to the lens of the 
line-of-sight to the source star and Vj_ is the relative transverse velocity of the 
encounter, both measured in the lens plane. The maximum optical amplification 
of the event is 



(4) 



The quantities A and T are easily derived from the event light curve, and 
constitute the physically significant observables for simple (i.e., point-source 
and point-lens) microlenses. 

Within a homogeneous set of observed microlenses, the overall statistical 
distribution of observables, /^,t {A, T), can be formally related to that of lens 
mass among the population from which they were drawn, /m (M), by a relation 
of the form 

fA,T {A, T) = J /^,r|M {A, T\M) /m (M) dM . (5) 

The notation is that fx {y) is the probability density function (pdf) of random 
variable x evaluated at a; = y; fx\z is the pdf of x conditional upon z. The 
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core of the model is the kernel function /^,r|M {A.,T\M), the probability that 
an object of mass M will yield a microlensing event characterized by {A, T) . 
This is essentially an accounting mechanism, one that adds up the probabilities 
of all possible combinations of lensing parameters {b,V±, Dl, Ds) that lead to 
the observed {A, T), given mass M. Its computation requires the provision of a 
kinematic model of the lens and source populations; i.e., a priori distributions 
of {b,V±, Dl, Ds) ■ Having computed the kernel, one numerically inverts the 
statistical model (eq. in terms of the observed /^,t to determine the desired 
mass function characterizing the population producing the microlenses. With 
kinematically well-characterized populations ~ e.g., those (presumably) of the 
Galactic halo and LMC ~ this procedure appears to be a reasonable one for 
estimating the mass function of the lensing population. 



2 STATISTICAL MODEL 

2.1 Derivation 

The detailed statistical relation can be derived as follows. Starting from the 
marginalization tautology 



/at (AT) 



fA,T,D^,Ds,M{A,T,DL,Ds,M) dDLdDsdM , (6) 



one changes variables within the integrand, from A and T to b — rs P (A) and 
Vi_ = rs/T: 

fA,r{A,T) = III h{rEP{A)) fv^M,,Ds,M(}Y^^^^^s,M) \3\dDLdDsdM 

(7) 

As shown in Appendix A, in any likely application the impact parameter b will 
be uniformly distributed independently of any other parameters, so that its 
distribution may be factored away from that of the other parameters and 
assumed to be a constant that can be absorbed into normalization of the model. 
The quantity |J| is the Jacobian of the variable transformation: 



db/dA db/dT 
dV^/dA dV^/dT 

Combining these results, 
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/^r (AT) = C- 
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\ fv^,D^,Ds,M ^i. Ds,M^ dDL dDs dM 

(9) 

where C is a constant and all of the amplitude information is contained within 
the likelihood function 
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which is independent of any other parameters in the problem. Equation ^ 
makes clear that A and M are effectively statistically independent, so the lensing 
amplitude carries no significant mass information. We can thus treat A as 
a nuisance parameter and integrate it from the model, leaving one in which 
the only observable information is that of timescale T. The integral of L^i 
actually diverges as ^ ^ 1, but that lower limit corresponds to unobservably 
faint microlenses: in any practical observing program there will be a minimum 
observable amplitude, however poorly defined, that serves to determine the value 
of the integral of equation ^ over A, and to partly convert the left-hand side 
of that equation to one of observed microlenses. To complete that conversion, 
we multiply the model by the probability E (T) that a microlens of timescale 
T will actually be observed; this may be estimated by the observer based on 
properties of the observing process, so the (non-normalized) model becomes one 
corrected for observational selection: 

fP {T)^^ J J J dD^dDsdM. 

(11) 

The superscript (O) on the observable fq- denotes a distribution of observed 
lenses only, while the (P) on the integrand refers to the distribution among 
the entire population of objects from which the observed lenses and sources are 
drawn. 

Put nearly into the standard form of equation [S] we have 

fP {'r)^C j Lr\M {r\M) /If > (M) dM , (12) 

where 

Lr\M {T\M) = J J {^^^^ D^, Ds\m) dD^dDs (13) 

is the likelihood of observing a microlens of timescale T arising from a lensing 
object of mass M (which enters via r^; cf. eq. [Q). This - equations 1121 and 
1131 - is the desired result, the deterministic connection between lens population 
mass function and observed timescale distribution (the normalizing factor C is 

to be computed a posteriori from normalization of (T))- 
2.2 A Corroborating Example 

As an example loosely based on observations of microlenses toward the LMC, 
we postulate a thick disk whose stellar spatial density varies along the line-of- 
sight to the LMC as p{D) cx exp (—£)//), where I = 8 kpc is the projected 
scale length. From Appendix C we then expect the linear distribution of stellar 
distances to be //^^ (Dl) oc Dj^exp{—D/l). We also postulate a Maxwellian 
velocity distribution for V± of 200 km s~^ rms, and put all the source stars at 
Ds = 55 kpc. Assuming statistical independence of Dl, Ds, V±, and M ; and a 
timescale observing efficiency of E (T) = T (200 — T) for < T < 200 days and 
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Figure 1: Relative likelihood ir|A/ ("^1-^) (eq- El? arbitrarily scaled) that a 
compact object of mass M will produce an observable microlensing event of 
timescale T, for a synthetic kinematic model and timescale observing efficiency 
(see text). 

otherwise; we use equation II 31 to compute the integral kernel shown in Figure 
1. With a presumed lensing population mass function uniform on (0,1) solar 
masses, this kernel yields (from eq. the timescale distribution of observable 
microlenses shown in Figure 2. Also shown in this figure are the corroborating 
results of one of many Monte Carlo simulations with sample lensing object 
kinematics drawn directly from the above distributions, without reference to 
the statistical model (the simulations employed 6max = 20 A.U. (Appendix A) 
and ^min = 1-3 for observable lenses, neither of which affects the results; and 
employs the same E (T) as used in the model). The excellent agreement between 
model predictions and an independent simulation suggests an accurate model, 
fully corrected for timescale- and amplitude-limited observational selection. 



The model may be numerically inverted to determine the underlying population 
mass function from an observed timescale distribution. The simplest inversion 
technique may be that of Richardson (1972) and Lucy (1974), in the form pro- 
posed by Heacox (1997, eq. [A4]). Suppressing the (P) and (O) superscripts 
for clarity, this iterative algorithm takes the form 



2.3 Solution 



(T, M) 



Lr\M {r\M) 



(14) 




/ Lr\M{A'm) fti^ (m) dmdt ' 



dt . 



(15) 



J J {t,m) f^Um) dm 
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Figure 2: Timescale distribution fr {T) for a hypothetical lensing population 
kinematic model and uniform mass function (see text). Solid curve: as predicted 
by the statistical model and computed from the kernel function of Figure 1 and 
adopted mass function. Histogram: as simulated by Monte Carlo sampling from 
the presumed kinematic distributions and mass function, without reference to 
the statistical model. 



One chooses an initial guess for flj' (M) (prudently chosen to be smooth on 

a wide range) and iterates until sensible convergence. If the observed is 
actually derived from the presumed kinematic model, the method will almost 
always converge to a reasonable approximation to the underlying mass function. 
Note, however, that the statistical model is a smooth, non-linear mapping from 
the space of mass functions into that of observed timescale functions, for which 
we do not expect that every timescale distribution will have a legitimate mass 
function pre-image, given the adopted kinematic model. If either the timescale 
function is inaccurate (duo, e.g., to sampling variations) or the adopted kine- 
matic model is incorrect, the Richardson-Lucy method may fail to converge 
when applied to this statistical model. 

3 APPLICATION TO GALACTIC HALO MI- 
CROLENSES 

3.1 Data 

The lensing population responsible for microlenses observed against the back- 
drop of the LMC probably represents the kinematically most homogeneous set 
of objects to be analyzed for mass function with the statistical model of the pre- 
vious section. The MACHO Project (Alcock et al. 2000) has detected between 
13 (Criteria A) and 17 (Criteria B) "simple" (point-source and point-mass) lens- 
ing events over a 5.7-year period. The timescale ((T)) distributions for these 
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Figure 3: Observed microlens timescale distribution toward the LMC. Solid, 

heavy curve: Alcock et al.'s (2000) high-confidence events ("Criteria A"); dashed 
curve, lower-confidence events ( "Criteria B" ) . Light curve: estimated timescale 
observability (a transcription of Alcock et al.'s (2002) 5-year Criteria A curve of 
their Figure 5). The symbols at the bottom denote the timcscalcs for individual 
Criteria B events (those of Criteria A are an approximate subset). 



observations are shown in Figure 3, together with their estimated observing 
efficiency E (T) (note that the timescale t of Alcock et al. (2000) is 2T) . 

Thetimescale frequency h.mctions as shown are Gaussian kernel estimators 
applied to Alcock et al's (2000) individual events (their Table 7), with a smooth- 
ing kernel standard deviation set as large as possible (7 days) consistent with 
the very low estimated probability of detecting events with timescales less than 
~ 2 days. The multi-modality of these distributions is quite possibly an artifact 
of small sample statistics, but the narrow range of observed timescales probably 
is not (the matter is discussed further in §4). 

3.2 Kinematic Model 

There appears to be no unambiguous kinematic model of the dark halo with 
which these data may be used in application of the statistical model. While 
there is some dynamical evidence for the dark halo mass distribution, there 
are few observational clues to the required velocity field even in the luminous 
halo. The root of the difficulty is that the velocity of interest here that on 
the plane of the sky - is unobservable beyond the solar vicinity and cannot be 
reasonably inferred from equilibrium dynamical models based on the collisionless 
Boltzrnann eq\.iation (Binney & Tremaine 1987). Beyond the solar vicinity, such 
velocities are largely transverse to the Galactic center and thus reflect orbital 
angular momentum as much as they do orbital energy. But angular momentum 
is essentially a free parameter in static equilibrium models, so that a whole 
sequence of models of different orbital angular momentum content may be flt to 
observed spatial densities and radial velocity distributions (Heggie & Hut 2003), 
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and there is thus no effective theoretical or observational constraint on halo 
transverse velocity distributions. Some evolutionary models of the luminous 
halo (e.g., Sommer-Larson et al. 1997) are consistent with reasonable levels 
of velocity anisotropy that, in conjunction with the observed radial velocity 
dispersion of distant halo stars, suggest a large (> 100 km s~^) transverse 
velocity dispersion throughout the halo. But it is not at all clear that such 
speculative dynamical models of the luminous halo can be confidently applied 
to the dark halo. 

One is thus reduced to examining the compatibility of presumed dark halo 

kinematic models with the observed microlcns timescale distribution. The ini- 
tial choice of kinematic model has here been made so that all parameters are 
statistically independent (the alternative will be considered shortly) and with 
the narrow range of observed timcscales in mind, so that the ranges of distances 
and velocities in the kinematic model have been kept as small as is reasonably 
consistent with a spheroidal, pressure-supported halo. Thus, the "NFW Law" 
of dark halo spatial mass density, p{R) oc (a + i?)~^ (Navarro, Frank & 
White 1997), is chosen as the basis of the MACHO distance distribution, where 
R is the galactocentric distance and a = 5 kpc is the core radius. The resulting 
solarcentric distance distribution is then (from Appendix C): 

fDADL)^Dl{R{DL)[a + RiDL)fY' , (16) 

1/2 

where R{D) = [D^ + i?| - 2Di?| cos 6lmc cos ?lmc] with Rq = 8 kpc 
presumed. LMC source star distances are taken to all be 55 kpc (the more 
realistic assumption of a depth to the LMC does not change the result). The 
transverse velocity V±_ of the lens relative to the source, as seen from our location, 
includes the motion of the Sun relative to the LMC, and of the non-rotating 
halo rest frame relative to the Sun. The first of these is computed from the 
observed LMC proper motions of /i„.(5 = 1.20,0.26 micro- arcscconds per year 
(Jones, Klemola & Lin 1994) and an LMC distance of 55 kpc, the second from 
a presumed rotation rate of the Sun relative to the Halo of 200 km s~^. The 
resulting mean velocity of the halo across the line-of-sight to the LMC is thus 
{V} = 215 km s"^. Relative to this mean frame we conservatively adopt a halo 
velocity dispersion on the plane of the sky of 120 km s~^ (Binney & Merrifield 
1998, Scheffler & Elsasscr 1988, Sommcr-Larscn ct al. 1997) and realize this with 
low angular momentum orbits characterized by a dark halo internal Gaussian 
random velocity distribution of mean 0. With these and the presumption of 
no significant LMC internal velocities, we compute (by the method of §3 of 
Appendix B) the distribution of V± shown in Figure 4. 

3.3 Mciss Function 

The integral kernel Lq-\M computed for these distributions (and assuming 

the E (T) estimated by Alcock et al. 2000 and illustrated in Figure 3) is similar 
in appearance to that of the example shown in Figure 1. In application to 
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Figure 4: Distributions of lens solarcentric distances Dl (left panel) and relative 
transverse velocities Vj_ (right panel) in the assumed dark halo kinematic model 
(see text). 

presumptive lens populations of a single mass (i.e., /m (M) as a delta function 
centered on the chosen mass), this model gives rise to the sample timescale 
distributions shown in Figure 5. 

It is apparent that no non-pathological lens mass function can adequately 
reproduce the observed timescale distribution, given this kinematic model (the 
Richardson-Lucy interactive inversion method of §2.3, applied to this model and 
Alcock et al.'s [2000] Criteria A observed timescale function, does not converge). 
The difficulty can be ascribed to the wide range of lens distances and velocities in 
the kinematic model, which lead to similarly wide ranges of observed timescales 
for given lens mass. 

The obvious alternative kinematic models do little, if anything, to rectify 
this discrepancy. Isothermal models, both with and without mass segregation, 
lead to different kernels Lq-\]^ (T|M) but similarly widely distributed timescale 
distributions irrespective of lens mass. A choice of different dark halo mass 
distribution (e.g., p (R) oc {a^ + i?^) or (a + i?)"^; Binney & Merrifield 1998, 
Alcock et al. 1996, 1997) will typically exacerbate the problem, as will the 
inclusion of LMC internal velocity dispersion in the computation of V±. The 
timescale distribution widths may only be adequately decreased in a spheroidal 
halo by (1) postulating a halo velocity field with a very small transverse velocity 
dispersion, probably less than ^ 50 km s^^, while continuing to ignore the LMC 
internal motions; and/or (2) greatly decreasing the spatial range encompassed 
by the lensing populations, as in a thin shell. Neither of these seems dynamically 
plausible. Additional caveats, based on the small sample size and poor numerical 
conditioning of the statistical model, are discussed in §4. For the moment, the 
microlens timescale distribution as observed appears to be inconsistent with a 
standard, spheroidal, pressure-supported model of the dark halo. 
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Figure 5: Microlcns time scale distributions (light curves) predicted for single- 
mass lensing populations from a spheroidal, pressure-supported, low angular 
momentum kinematic model of the halo (see text). The distributions arc labeled 
with the stellar mass (solar masses) that produced them. The heavy curve is 
the adopted observed timescale distribution (Criteria A of Alcock et al. 2000). 



3.4 Thick Disk Model 

As one possible alternative to the dark halo microlensing model we consider 
a rotating, low-dispersion, thick disk as the site of microlensing objects. To 
be consistent with the roughly constant rotation velocity in the outer disk, we 
adopt a spatial mass density that varies with galactocentric distance r in the 
disk, and height z above the disk, of p (r, h) oc exp {—z/h). Following the 
argument of Appendix C, the distribution of solarcentric distances will be 

for. {DL)cxDlr{DLr^e^p 

where r (D) = [D^ cos^ &lmc + R% — 2DRq cos Blmc cos Ilmc] and we again 
choose Rq = 8 kpc. To keep the timescale dispersions small we choose h = 2 
kpc in a co-rotating thick disk of 30 km velocity dispersion. Both of these 
are at the low end of what would be expected, both observationally and dynam- 
ically; but they seem possible. We further assume the statistical independence 
of all kinematic parameters, a constant source distance of 55 kpc, and an in- 
ternal velocity dispersion of 30 km s ~^ for the LMC (but no rotation). Then 
the integral kernel Lt|m C^I-^) corresponding to this kinematic model produces 
the single-mass timescale functions of Figiirc 6, from which it is apparent that 
a distribution of lens masses encompassing roughly 0.1 to 0.7 solar masses may 
suffice to explain the MACHO observed timescale distribution. In this case the 
Richardson-Lucy inversion procedure of §2.3 converges to the results shown in 
Figure 7. For this postulated thick disk there is no evidence for lensing objects 
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Figure 6: As with Figure 5, but for a rotating, low velocity dispersion, thick 
disk population (see text). The thin curves are labeled with the lens mass (solar 
masses); the thick curve is the MACHO observed timescale distribution (Alcock 
et al. 2000). 
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Figure 7: Thick disk mass functions approximate solutions. Left panel: mass 
functions (arbitrarily scaled) inferred from the Richardson-Lucy algorithm ap- 
plied to the statistical model derived from an adopted thick disk kinematic model 
(see text) and observed timescale function (solid curve of Figure 3). Solid curve: 
200 interactions beginning from a uniform mass distribution. Dashed curve: 100 
interactions. Right panel: the timescale distributions (light curves) correspond- 
ing to the mass functions in the left panel, as computed from the statistical 
model (eq. ^2] and I13|l . The heavy curve is the observed timescale distribu- 
tion. 
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exceeding one solar mass, and some evidence that a significant fraction of the 
lenses have low- or sub-stellar masses. 

4 SUMMARY AND CONCLUSIONS 

The statistical model of equations El and El is the complete deterministic 
connection between the lens population mass function and the distribution of 
timescales among observed microlenses. It is fully corrected for observational 
selection based on timescale and maximum amplitude observability, providing 
these are mutually independent (which will usually be the case) . The correction 
for selection on the basis of amplitude arises from a basic model (eq. j^l) that 
includes the amplitude as an observable, from which it is removed by integra- 
tion over, in effect, observable amplitudes only. The result (easily proven by 
construction of a basic model for timescales only) is the inclusion of an addi- 
tional factor of in the integrand: the observability of a microlensing event, 
in terms only of amplitude, is proportional to the Einstein radius of the lens, 
hence (approximately) to M^/^. 

Appendix A provides a formal derivation of event impact parameter distri- 
butions: for all practical purposes these are uniform on a range exceeding that 
likely to be encountered in amplitude-limited observations. These results lead 
to a formal distribution of event amplitudes (eq. ^01): which is independent 
of kinematic model parameters and of lens mass; so that microlens mass infor- 
mation must come only from event timescales. These are both familiar results, 
here rigorously derived. 

The statistical model is limited in applicability by its relatively poor numer- 
ical conditioning, reflective of the small contribution of lens mass to the event 
timescale. The formal concept of condition number is not strictly applicable 
to such non-linear models, but the sensitivity of the observed timescale distri- 
bution to the mass function can be crudely estimated as the change in mean 
timescale (as computed from the model) for a given change in scale of a mass 
delta function. For the example of §2.2, this sensitivity is about A (T) « 7 
days per AM = 0.1 solar masses over the range of (0.1, 1) solar masses. While 
this number depends upon the kinematic model adopted, it is probably typical 
of most applications to widely distributed kinematics and suggests a relatively 
poorly conditioned model from which mass function information will be difficult 
to extract. 

The poor numerical conditioning of the model is reflected in its smoothing 
effect in translating mass functions into timescale distributions. This is appar- 
ent from the thick disk model result of §3.4, where nearly delta- function mass 
distributions serve to reproduce the smooth, bimodal observed timescale dis- 
tribution. It is for this reason that no realistic mass function can reproduce 
the substantial timescale structure in Alcock et al.'s (2000) Criteria B sample 
(dashed curve in Figure 5): the isolated event near 150 days, if real, is a true 
outlier, perhaps arising from a separate lens population. 

For microlensing toward the LMC, the conditioning difhculty is exacerbated 
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by uncertainties in the observed timescale distribution due to small sample 
statistics. Re-sampling of the presumed timescale distribution (Criteria A of 
Alcock et al. 2000; heavy curve in Figure 5) demonstrates that a wide range of 
timescale distributions is likely to arise from such a small sample. A significant 
fraction of 13-object distributions drawn from, e.g., a Gaussian of mean 50 days 
and standard deviation 15 days (roughly those of Criteria A events) show bi- 
modality similar to that of the Criteria A distribution, so that the observed 
bi-modality may well be a spurious consequence of small-number statistics. 

It thus seems unlikely that the current microlensing observations toward 
the LMC could meaningfully constrain the (presumably) underlying dark halo 
mass function. Even such summary statistics as the mean mass are in some 
doubt: since M oc for given event timescale, and the mean of V± depends 
as miich upon the poorly known LMC proper motions as it does on the solar 
galactocentric velocity, the canonical value of (M) 0.5 solar mass (e.g., Alcock 
et al. 2000) may well be incorrect by a substantial factor. 

But the narrow range of observed timcscalcs appears to be real: re-sampling 
experiments from this and similar distributions suggest little likelihood that the 
true distribution is sufficiently broader than that observed as to be compati- 
ble with kinematic distributions as broadly distributed as those expected of a 
spheroidal, pressure-supported population of modest mean orbital angular mo- 
mentum. It seems likely that the population producing the microlenses observed 
toward the LMC must come from something other than a "standard" halo, ci- 
ther dark or luminous. The results of §3.4 suggest that a thick disk origin is one 
possibility; no doubt there are others. 
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A IMPACT PARAMETER DISTRIBUTION 



Wc take as the most likely geometric model one in which, at the commence- 
ment of the survey, the lensing candidates are uniformly distributed on the 
plane of the sky in the vicinity of the source star being observed. The geometry 
of the encounter is described entirely by the original distance between lensing 
object and source star, r, and the direction ^ of motion relative to the radius 
vector from star to lens; both measured in the lens plane. The impact param- 
eter corresponding to (r, ^) will then be 6 = r sin ^ ; r will be distributed as 
fr{r) = 2r/rj^jj^ where Tmax is the limiting radius of the search area in the 
lens plane; and, by virtue of the uniform distribution of lens locations, we can 
safely presume ^ to be uniformly distributed on (0,7r/2). By a useful theorem 
of cumulative distributions (e.g., §7.26 of Kendall, Stuart and Ord 1987), the 
pdf of b is given by 



fr{r)f^ix) d^dr 



(18) 



where £7 = {('":'?) |?'sin^ < 6} is the set of all lens coordinates for which the 
impact parameter is less than b. Thus: 



h{b) 



d 

'db 



b /'7r/2 r^msLX /'arcsin(6/r) 



b JO 
1/2 



4r 



d^dr, (19) 



(20) 



for < < r-jjiax and otherwise. This formula has been verified by Monte 
Carlo simulations based on the above geometric model. In most applications 
the largest impact parameter corresponding to observable amplitudes will be 
smaller than rmax by several orders of magnitude (about 6—10 for observable 
halo microlenses) and fb {b) may safely be taken to be a constant independent 
of all other parameters. The above result, with the assumption b <C f'maxj leads 
to the amplitude likelihood of equation^] From that relation the distribution 
of scaled impact parameter /3 is easily seen to be 



fpm^LA (Am 



dA 



dp 



= 1 



(21) 



The distributions of observable microlens impact parameters, scaled or not, are 
thus expected to be uniform and independent of distributions of lens distances, 
velocities, or masses. 



B VELOCITY DISTRIBUTIONS 

The velocity Vj_ is the relative velocity of lens and source star, projected onto 
the plane of the sky, as seen from our location. It can usually be computed as the 
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magnitude of the vector sum of two or more projections of velocity components 
onto the plane of the sky. For many applications, the following results will 
suffice to compute the distribution of this quantity, from a kinematic model of 
the three-dimensional velocity distributions of lenses and source stars. 

B.l Projection onto Plane of Sky 

An isotropic space velocity V will have a magnitude projected onto the plane 
of the sky, Vpos, given by the simple projection relation encountered in, e.g., 
analysis of binary star orbits (Chandresakhar and Miinch 1950): 

oo 

(Vpos) = Vpos I f""^^^ dV . (22) 



B.2 Combination of Random Velocities 

Two statistically independent, mutually orthogonal, random velocity compo- 
nents {Vx, Vy) on the plane of the sky; have a combined mag nitude V = [V^ + V^] 
at a position angle 9 — arccos (14,/^), whose joint distribution is given by the 
simple variable transformation: 

fvfi {V, e) = V fv^ {V sin 9) fv, {V cos 0) (23) 

(V is the Jacobian). The marginal distribution of velocity magnitude alone is 
given by the integral of this expression over or, changing variables to Vx = 
y sin 61, 

}fvAV.).fy(VV^^) 
fviV)^vJ ^^==-^dVx. (24) 

B.3 Combination of Random & Streaming Velocities 

The distribution of magnitudes of the combination of a random velocity Vr with 
a streaming velocity Vs (e.g., that arising from the LMC proper motions), both 
on the plane of the sky, may be computed as follows. Let fva.o {Vr,0) be the 
joint distribution of random velocity magnitudes and directions (as measured 
from the streaming velocity direction), as computed by equation 1231 or. e.g., 
as a Gaussian velocity whose mean and/or standard deviation is a function of 
direction. Then, following the same procedure as with eauation l24l the resulting 
velocity magnitudes on the plane of the sky will be distributed as 

f n/ ^ 1/ J fvn.e{VR{V^,<j>),0{V^,cj,)) ^^ 



1/2 
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where 

Vr{V^,<P) = [Vi + Vl-2VsV^coscpY^^ 

V± sin (j) 



sin (0(^^,0)) 



and V±/Vii {V±, (f) is the Jacobian of the variable transformation from (Vr, 0) 
to (Vj_, (j)). If the random component is isotropic, this reduces to the relatively 
simple expression 



fvAV^) = ''-^ I ^^^^ dVn ■ (26) 

^ ^J{2VsV^f-{V| + Vl-V^f 



|Vs-Vj.| 



An example may be seen in Figure 4. This rather unusual function has been 
verified by Monte Carlo simulations independent of the statistical reasoning 
underlying this derivation. If Vg — * this expression reduces to fv^^ = Jvr , as 
it must. 



C DISTANCE DISTRIBUTION 

If the spatial number density of stars varies with galactocentric distance R as 
p{R), and the distribution is radially symmetric about the point from which 
R is measured, the linear density (number per unit change in R) is easily seen 
to be fR{R) cx R^p{R). But viewed from off-center, or if the spatial distri- 
bution is not radially symmetric, the matter is complicated by the dependence 
of the density upon the viewing direction. In the case of microlenscs viewed 
against the LMC, the field of view is considerable and it is not obvious how 
the spatial number density will translate into a linear density in solarcentric 
distance D when viewed from our off-center position. Formally the solution is 
straightforward: by the usual change of variables, 



fD,b,i {D,b,l)<xR^p{R) 



d{RA,0) 



d{D,b,i) 



(27) 



1 /n 

where R = [Rq + — 2RqD cosbcos and {cl>,e) are the galactocentric 
angular coordinates for which we presume connecting relations with the solar- 
centric position (e.g, sin^ = D sinb/ R{D,b,l), etc.). The average value for the 
field of view is just the integral of this density over the corresponding spreads 
in b and /. The Jacobian of this relation is difficult to represent analytically, 
but is numerically trivial. As the field of view shrinks the spatial density be- 
comes more uniform over that field until, in the limit of a small field of view, 
the solarcentric distance pdf becomes very nearly 

fD,b,i {D, b,l)ocD^-p [R {D, b, I)) , (28) 
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where (6, 1) refer to the center of the field. In numerical experiments for the 
LMC direction, with a ~ 10° wide field and an NFW spatial density profile (cf. 
§3.2), this approximation agrees with the more exact result (eq. [23) to within 

better than 1%, when integrated over the field. 
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